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Abstract 

Existing calculations of heavy quark production in charged-current and neutral current lepton-hadron scat- 
tering are formulated differently because of the artificial distinction of "light" and "heavy" quarks made in 
the traditional approach. A proper QCD formalism valid for a wide kinematic range from near threshold to 
energies much higher then the quark mass should treat these processes in a uniform way. We formulate a 
unified approach to both types of leptoproduction processes based on the conventional factorization theorem. 
In this paper, we present the general framework with complete kinematics appropriate for arbitrary masses, 
emphasizing the simplifications provided by the helicity formalism. We illustrate this approach with an 
explicit calculation of the leading order contribution to the quark structure functions with general masses. 
This provides the basis for a complete QCD analysis of charged current and neutral current leptoproduction 
of charm and bottom quarks to be presented in subsequent papers. 
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1 Introduction 



Total inclusive lepton-hadron deep inelastic scattering has been the keystone of the quark-parton 
picture and the QCD-based Parton Model. As the global QCD analysis of high energy interactions 
becomes more precise, other processes be^ii^tixplay an increasingly important role in determining the 



e, other processes be^iiLta.pia; 
inside the nucleon.B B31 



parton distributions inside the nucleon.lir cr B ul For instance, semi-inclusive charm production in 
charged-current and neutral-current interactions in lepton-hadron scattering serves as a unique probe 
of the strange quark and charmed quark content of the nucleon.B u LP In general, the production of 
heavy flavors in lepton-hadron and hadron-ha 



QCD study and for searches for new physics! 



™ rn l lirlpi-s is a yprv inr 



j^tant tool for quantitative 

Traditional analysis of massive quark production in DIS uses the simple light flavor parton model 
formulas (based on tree-level fpxwardXlpmpton scattering off the quark) with a "charm threshold" 



formulas ( based on tree- level ipxw.am.ijpmpton scattering on tne quark j witn a charm threshold 
or "slow-rescaling" correction JtZI E3 This prescription is still widely used-.in, current literature, 
particularly for dimuon production in neutrino charged current scattering.B o LJ; however, the ap- 
plicable range of this approach is very limited - for the neutral current case by the mass of the initial 
state quark; and for both cases, by the numerically important next order gluon contribution. KiJ In 
most neutral-current charm production calculations and recent HERA studies of heavy flavor produc- 
tion, a contrasting view has been provalent: one forsakesJaading-order quark scattering mechanism 
and concentrates on the 0{a s ) "gluon-fusion" processesMir c3 Whereas this latter approach is ap- 
propriate when theJiard scattering scale of the process, say Q, is of the same order of magnitude as 
the quark mass m , it is a poor approximation at high energies. In fact, when m/Q is small, these 
"gluon-fusion" diagrams contain large logarithms, i.e. factors of the form a™ log™ (m/Q), which viti- 
ates the perturbation series as a good approximation. These large logarithms need to be resummed, 
which then yield quark-scattering contributions with properly evolved parton distribution for the 
not-so-heavy massive quark. 

A consistent QCD analysis of this problem requires a renormalization scheme which contains the 
two conventional approaches as limiting cases — in their respective region of validity — and provides 
a smooth transition in the intermediate region where neither approximation is accurate. Such a 
scheme, motivated by the Collins- Wilczek-Zee cff renormalization procedure, was proposed some 
time ago in the context of Higgs production, resulting in a satisfactory theory valid from threshold 
to assymptotic energies. E^T This approach also provides a natural framework for heavy quark 
production. It is particularly simple to implement in leptoproduction-production processes, and has 
been applied to charm production in DIS in a previous short report .EU 

The current paper is the first of a series which will give a detailed formulation of this problem. In 
systematically developing a consistent formulation of heavy- flavor production in DIS, one finds that 
conventional calculations, even at the leading order level, make implicit approximations inherited 
from the zero-mass parton model-such as the Callan-Gross relation and the choice of the scaling 
variable-which are not always valid in the presence of masses. In order to make a fresh start on 
a consistent theory including non-zero-mass partons, this first paper is devoted to a self-contained 
development of the general formalism of deeply inelastic scattering in the presence of masses which 
is valid for both charged and neutral current interactions. Much of this is kinematical in nature. In 
considering charm production in existing fixed target neutrino experiments, an important practical 
consideration is that the target nucleon mass is comparable to the charm quark mass, and both are 
non-negligible compared to the average energy scale Q of the process. Thus, for consistency, target 
mass effects should also be incorporated precisely. To this end, we present a helicity formalism (along 
with the conventional tensor approach) to develop the general framework. It will become clear that 
whereas the conventional tensor method becomes quite complicated when both target mass and 
quark mass effects are properly incorporated, the helicity formalism retains the same simplicity 
throughout — due to its group-theory origin, and to a key feature of the QCD Parton Model. To 
make the general formalism concrete, we shall apply this helicity approach to a complete leading 
order calculation of heavy flavor production in charged current DIS, and then compare with the 
conventional tensor calculations. Numerical studies will show that the complete calculation (with 
all masses retained) leads to significant differences in the calculated cross-sections in certain regions 
of phase space. In the text of this paper, we shall emphasize the key elements of these developments. 
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Tab. 1: The gauge couplings of the vector bosons according to the Standard Model. 



Most technical details are relegatcd-ip the appendices. 

The second paper of this series shall be focused on the consistent QCD formulation of heavy 
quark production in the context of order a s calculation of this process, using the general kinematical 
formalism developed here. The emphasis will be on the formulation of a consistent renormalization 
and factorization scheme to reconcile the quark-scattering and the gluon-fusion mechanisms. The 
QCD framework developed there applies to all heavy quark processes, including hadroproduction. 
In subsequent papers, we shall study the phenomenological consequences of these calculations on 
the analysis of existing dimuon data from fixed target experiments, and on predictions of charm and 
bottom production at HERA. 

2 Scattering Amplitudes 

We consider a general lepton-hadron scattering process^] 

£i(£i) + N(P)^£ 2 (£ 2 )+X(P x ) (1) 

as depicted in Fig. [j] where the exchanged vector boson (7, W, or Z) will be labelled by B and its 
momentum by q. 




Fig. 1: The general lepton-hadron scattering process: N(P) + £1 — > X(Px) + £1 via the exchange 
of a vector boson, B{q). The lepton momenta are £i while the initial and final hadronic 
momenta are P and Px, respectively. 



The lepton-boson and quark-boson couplings are specified by the following generic expression for 
the effective fermion-boson term in the electro- weak Lagrangian: 



EW 



-9b 



(2) 



where a summation over B is implied. The gauge coupling constant gs for the vector boson field 
Vb depends on B and their values as prescribed by the Standard Model are given in Table |. 

+ X' where X' is unobserved. For the 



2 In the production of a heavy quark Q, the final state is given by X 
purposes of the present discussion, we shall not single out Q from X. 
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Both the hadronic and fcrmionic current operators are defined by 



4 W = t.f( x )ln (9v ~9Al b )V>/( 



x 



— (31 
= ^ f {x) ltl [g R (l + 7 5 ) + 9l{1 - 7 5 

where ipf denotes a generic fermion field, and the vector and axial vector couplings gv,A are related to 
their chiral counterparts by gL,R by gv,A — gL± 9r- The values of those fermion coupling constants, 
according to the Standard Model, are given in Table ||; however, we will keep them general in our 
considerations. 

The scattering amplitude for the process of Eq. ([!]) — with particle momenta as shown in Fig. [I] — is 
given by 

M = r ^ q) Wrm/ M (4) 

where q = - i 2 , I = l x + h, Q 2 = -q 2 > 0, and G% = g» v - q^q v /M%. The lepton current 
matrix element is given by 

= (h\r\h)^u(h)r[9B,(i + i 5 )+gL(i-i 5 )Hh) (5) 

The hadron current matrix element is kept in the general form: J*(P,q) = (Px \ \ P) ■ For simplicity, 



we have suppressed the polarization indices for all external particles in Eq. (|4|). Furthermore, the 
term can be replaced by g^ in actual applications since the term proportional to q^q v (when 
contracted with the lepton current matrix element) yields terms proportional to mj/Q 2 which are 
negligible at high energies. 

An alternative expression to the above familiar formulation of the scattering amplitude which 
emphasizes the helicity of the exchanged vector boson is given by:[^7], 

m = j; n (Q 2 ,p- q ) g \flZi l nQ2) (6) 

where n and m are helicity indices for the vector boson, j n {Q 2 ) and J^Q 2 , Pq) are the scalar helicity 
amplitudes for the two vertices shown in Fig. []J and d}(ift) is a spin 1 "rotation" matrix sp ecif ying 
the relative orientation of the two vertices. The derivation of this formula can be found in ¥3 c3 ; 
the precise definition of the rotation angle^J jb is given in Appendix EL (See also Appendix [b| for 
details). We note that the structure of Eq. (||) is quite similar to Eq. (Q) above. The advantages of 
using the helicity formulation in the QCD analysis of heavy quark production will be discussed in 
Section |. 

3 Cross-section Formulas and Hadron Structure Functions 

The cross-section formula for this process is (c/. Appendix |A|), 

d(T = 9Ar Gl ? 2 A^ 4 *Q 2 L " K dT (7) 

where Gi = g 2 B / (Q 2 + M^ ) is a short-hand for the boson coupling and propagator. The two indices 
B\ and Bi denoting the species of the exchanged vector bosons are implicitly summed over and kept 
distinct to accommodate the possibility of i~Z interference, and dT is the phase space of the final 
state lepton. The factor AttQ 2 is from the normalization of L and W. In the above expression we 
have introduced the dimensionless lepton and hadron tensors given by^] 

L ^ = iE<4b#2>^|f|fi) (8) 

spin 

w ^ = ^E{^Y$\P + q-Px)(P\J»\Px)(Px\Jt\P) (9) 

spin 



3 For space-like q, tp is actually a hyperbolic angle specifying a Lorentz boost. 

4 Historically, the definition of — and thus the definitions of Wi in Eq. ( px| ) — contains an extra factor of M, 
the target mass. In view of scaling considerations, it is more natural to use the dimensionless definition. Also note 
that sums and integrals over all the unobserved hadronic final states X arc implied in Eq. (H). 
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Tab. 2: The gauge couplings of the vector bosons according to the Standard Model. Vij represents 
the CKM flavor mixing, if relevant, and Qi is the fermion charge in units of |e|. 



The explicit expression for L M „ with general coupling constants is given in Appendix [b| As is well 
known, the hadron tensor W^ v can be expanded in terms of a set of six independent basis tensors^] 



P^Pu 



<tqu 
' M 2 



Wa 



M 2 
P»q v 



Wo 



2M 2 



1 v -w, 



■w 3 + 

P»q v 



(10) 



2M 2 



W 6 



where M is the target mass and e P9AIly = e"^^ P a q/3- The scalar coefficient functions {Wi} are the 
invariant hadron structure functions for this process. 

By substituting the lepton and hadron tensors in Eq. (0) and partially integrating over the phase 
space of the final state lepton one obtains, in the limit of negligible lepton masses, the well-known 
cross section formula, generalized to arbitrary couplings, 



da 



dE% d cos 6 



2E 2 Gi G2 
ttM ng 



2Wx sin 2 



W 2 cos 2 



±9 2 . 



Ei + E-2 
M 



Wa sin 2 ■ 



(11) 



where the ± sign for the W3 term refers to the case of lepton/anti-lepton scattering, respectively. 
Here, E\ and E2 are the energies of the initial and final state leptons respectively in the laboratory 
frame, 9 is the scattering angle of the lepton in the same frame, and ng is the number of polarization 
states of the incoming lepton. To simplify the expression, we define g\^ = g\ t ± g\g, where gLt 
and gn 1 refer to the chiral couplings of the vector boson to the leptons[] 

It is worth noting that the hadron structure functions {W4, W5, We} do not appear on the right- 
hand side because they are multiplied by factors of lepton mass from the lepton vertex, not because 
they are intrinsically small compared to the familiar {Wi, W2, W3}. This will become relevant when 
we discuss the calculation of hard scattering cross-sections involving heavy quarks. 

5 In some papers, the tensor associated with W\ is chosen to be the gauge invariant form (— g„ + q^q^/q 2 ), and 
that associated with W2 is obtained with the substitution P M — » P v ((jv — q^qv/q 2 )', these changes (convenient for 
conserved currents) will modify the definitions of W4, W5 and Wa only. 

6 The lepton chiral couplings appear explicitly because Ljl has been evaluated. The corresponding hadron chiral 
couplings reside in the {Wi} invariant structure functions. 
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It is by now customary to introduce the scaling structure functions Fi given by 





= Wx 
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M 




V 

= — w 3 
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(12) 



in terms of which the expression for the differential cross section may be rewritten as 



da 
dxdy 



2ME 1 GxG 2 
7T ni 



2 

9+£ 



xFi y 2 + F 2 



(i-y) 



(Mxy\ 
\2E 1 J 



±gi i [xF 3 y(l-y/2)} 



(13) 



In the alternative helicity formalism, the expression for the cross section is given by 
da yQ 2 GiG 2 f 2 \(F + + F_ 



dxdy 



2- 



ii t 



2 

9+t 



-(1 + cosh ip) + Fq sinh ip 



T9±A(F+-F-)coshiP] 



(14) 

where ip is the hyperbolic rotation angle of Eq. (^|), and we have introduced the helicity structure 
functions {F\, A = +,0, — } which correspond to the physical forward Compton scattering helicity 
amplitudes 



F x = e**(P,q) W v (P,q) e v x (P,q) 



(no sum over A) 



(15) 



with right-handed (+), longitudinal (0), and left-handed (-) vector bosons respectively.^ We note that 
the first term on the right hand side involves the transverse structure function Ft — {F+ + FJ)/2, 
whereas the third term is the parity- violating term with F + — F- proportional to F% in Eq. (|l3|). 
Eq. (14) should be familiar, as it is analogous to the corresponding well-known formulae for time- like 
vector boson production processes — Drell-Yan pairs and W-, Z-production — where the hyperbolic 
angle ip is replaced by the center-of-mass angle 6 for the final state lepton pair. 

The helicity structure functions as defined above are naturally scaling functions. In addition, 
their direct physical interpretation leads to simple properties in the QCD parton model framework, 
as we shall see in the next section. Note that Eq. (|lj) does not show any explicit target mass 
dependence; all complications arising from the non- vanishing mass are contained in the definition of 
the rotation angle ip through kinematics. This simplicity is a consequence of the underlying group- 
theoretical approach to the factorized structure of Fig. [j]. The precise relations between the helicity 
structure functions and the invariant structure functions are found (c/., Appendix |b|) to be: 



F+ = 
F_ = 

= - 
that Q 2 /v 2 



Fi 
Fi 
Fi 



1 



14 
1 + 3 



(16) 



(£) F2 



and we obtain the approximation: F± ~ Fi =F F%/2 and 



We see in the limit M 
F ~ -F 1 +F 2 /2x. 

To leading order in the electroweak coupling, Eq. (|l^), Eq. ( p^ ) and Eq. ( pT| ) are completely 
general, assuming only Lorentz kinematics and small lepton masses. In particular, all results up 
to this point are independent of strong interaction dynamics. Aside from Eq. (|l4|), they are well- 
established formulae explicitly generalized to include arbitrary couplings. 



4 The QCD Factorization Formulas 

Perturbative QCD allows one to relate the measurable hadron structure functions {Fi} to the cor- 
responding quantities involving elementary particles — the partons — which can be calculated in per- 
turbation theory. This section states the basic QCD "factorization theorem" as it applies to deeply 

7 The choice of these labels — over the more obvious R, L, etc. — is constrained by the conflict between the Left- 
handed and Longitudinal designations. For mg = 0, we can ignore F\ = {F qq , F q o, Fo q }, cf., Appendix H. 



6 



inelastic scattering processes and points out some important unfamiliar features in the presence of 
non-zero masses, especially when the initial state parton is a heavy quark. 



4.1 Factorization of Tensor Amplitudes 

The factorization theorem states that, in the Bjorken limit, the dominant contributions to the 
hadronic tensor structure function has the factorized form of Fig. || with on-shell, collinear partons: 




Fig. 2: Pictorial representation of the factorization theorem for the hadron structure functions for 
inclusive deeply inelastic scattering. The process on the left is N(P) + B(q) — > X(Px), and 
the factorized process on the right is N(P) — > a(k±) (represented by the parton distribution 
function, /jy) with the successive hard scattering interaction a(fei) + B(q) (represented by 
a;" ). The vertical lines indicate an inclusive sum over the final states, X(Px)- 

In Eq. (|l7|), the label 'a' is summed over all parton species. The convolution integral variable £ is 
the momentum fraction carried by the parton with respect to the hadron defined in terms of the ratio 
of light-cone momentum components £ = / P + . The universal parton distribution functions 
are scalars; scattering of the vector boson takes place with the partons via the hard-scattering factor 
o^f" which can be aptly called the parton structure junction tensor since it is entirely analogous 
to the hadron structure function tensor W^ N by substituting the hadron target W by the parton 
target 'a! Note, the tensor structure of W/H is completely determined by that of oj^. These features 
should be obvious by inspection of Fig. 0. Strictly speaking, the factorization theorem is established 
in this simple form only for certain specifically defined asymptotic regimes. We shall treat Eq. ( |l7j ) 
as an ansatz and apply it in such a way that our results reduce to the known correct expressions in 
the limits A <C m,2 — Q on the one hand, and A < iri2 <C Q on the other. 

The presence of heavy quarks among the initial and final state partons in lu^ has some important 
consequences. The most immediate one is that the range of integration in Eq. (|l7]) will depend on 
the masses of the heavy quark as a simple consequence of the kinematics of the hard scattering. 
In leading order QCD, where the integration range reduces to a single point, this naturally gives 



rise to a generalized 
simple parton model 



ow-rescaling" variable which was originally proposed in the context of the 
(Cf., Appendix |a|.) In addition, the tensor structure of the perturbatively 
calculable lo^ is clearly different from that of the naive parton model, even in leading order QCD! 
For example, the well-known Callan-Gross relation simply does not hold in the presence of heavy 
quark mass. A proper treatment of heavy quark production must use the correct hard-scattering 
amplitude uj^ (calculated to the appropriate order, including quark masses) in conjunction with 
choosing the proper variable. A "slow-rescaling prescription" of a simple variable substition is not 
sufficient, cf., Section 0. 

In order to apply the factorization theorem to measurable quantities properly, we must re-express 
Eq. ( |l7j ) in terms of the independent invariant structure functions {Wi} or the helicity structure 
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functions {F\} in a precise way. Theoretical calculations of the parton-level hard amplitudes on 
the right-hand side of the equation usually yield the (parton) invariant or hclicity amplitudes, not 
the tensor uj 111 ' itself. In the presence of target and heavy quark masses, we will find that the 
relations between the invariant structure functions at the hadron and the parton levels are far from 
being simple, as usually assumed in existing literature. In contrast, the connection between the 
corresponding helicity structure functions are completely transparent. 



4.2 Invariant Structure Functions: 

The parton-level invariant amplitudes uii are defined in analogy to Eq. ([l^), as follows:^ 



q*q v k^q„ + q^k lu k^q„-q^k lv 

+ W + 2Q5 W6 (18) 



where k% is the momentum of the incident parton. Substituting Eq. (|lq ) in Eq. ( |17[ ) and comparing 
uj^ u with W^ v (Eq. ([To|)), we see that the relations between invariant structure functions at the 
hadron and the parton levels depend on the relation between fc^ 1 and P M . Whereas the two momenta 
are proportional in the zero mass limit, this relation becomes non-trivial in the presence of either 
target mass or parton mass, (c/., Appendix [A]) . Since the vectors P, ki and q are collinear, we can 
parametrize k\ as 

K = CpP» + C9" (19) 

In the zero mass limit, £p — ► £ and £ 9 —* 0. In general, the coefficients (£p, ( q ) are rather complicated 
functions of the masses and the convolution variable £, (c/., Eq. (|99|)). Thus, the relations between 
the Wi and the Ui are also rather complicated. Relevant formulas which relate Wi to w, are given 
in Appendix IB. 



4.3 Helicity Structure Functions: 

In sharp contrast to the above, the factorization theorem assumes a simple form when expressed in 
terms of the helicity basis. To see this, let us define the parton helicity structure functions u>\, in 
analogy to Eq. (|l5|), by: 

"A = el*{ki,q)bJ lx v e' / x {ki,q) (no sum over A) (20) 

In order to relate these to the hadron helicity structure functions F\, Eq. (|l^), it appears that one 
needs to re-express the vector-boson polarization vectors {e^ (k± , q)} (defined using k\ as the reference 
momentum) in terms of {e^(P, q)} (defined using P as the reference momentum). The enormous 
simplification of the hclicity approach follows from the fact that the two sets of polarization vectors 
are in fact identical even in the presence of masses, hence no transformation is needed! The reason 
for this is that for a given vector-boson momentum q, the reference momentum is used only to 
specify the direction of the polarization axis; the two seemingly different reference momenta k\ and 
P actually specify the same set of polarization vectors because they are collinear in the QCD Parton 
framework. Thus, we arrive at the straightforward formula: 

F* N (<l,P,-)=Y,fN®"x' a (21) 

a 

This suggests that to explore the consequences of perturbative QCD on heavy quark production 
(as well as on all other processes), it is advantagous to perform the calculation in the helicity basis. 

8 In order to render the uii dimensionless, we use the natural variable Q rather than any parton mass in scaling 
the tensors so that the invariant structure functions have well defined limits as m/Q — » 0. (Note, if the hadronic 
structure functions were originally denned this way, rather than using the target mass M as the scale factor, {W^} 
would be naturally "scaling!") 
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The simple formula Eq. (^T]), together with Eq. (|l4|), relate the calculation of hard scattering ampli- 
tudes directly to measurable cross-sections without any approximations or complications. Besides, 
since the parton-level helicity amplitudes have simple symmetry and structure, due to the basic 
chiral couplings of the theory, the results of this approach are often the most physical and compact 
to begin with. 



5 Leading Order QCD Calculation of Heavy Flavor Production 

To illustrate the use of the general formalism developed above, we apply it to the calculation of 
heavy quark production in leading order QCD. Existing applications of heavy quark production in 
DIS mostly concern charm production in charged current interactions at fixed-target energies. Since 
the charm mass is comparible to the target mass for existing neutrino experiments, and neither 
is negligible compared to the energy scale Q, it is reasonable to retain the target mass effects in 
order to be self consistent. Numerical comparisons of the complete calculation (with full target mass 
dependence) to the conventional one show that the difference can be significant in certain regions of 
the phase space. 




Fig. 3: Leading order hard-scattering amplitude for heavy quark production. 

The leading order diagram that contributes to uj\ is shown in Fig. ^ and its contribution, including 
all masses and arbitrary couplings, is calculated explicitly in Appendix |c[ We consider charm 
production in charged current neutrino scattering. Since, the W-exchange process involves only 
left-handed chiral couplings, (c/., Table |^). The parton helicity structure functions for scattering 
from a strange quark are given by 

12 



w ± = 91 a 



2 Q J + mj + mjTA s ( j 1 

T. „ T I — 1 

A \x 



2 (mj - mXflQ 2 + m|+mf / j . 
wo = 9 La x S \ 1 ( 22 ) 



A \ X 

where g 2 is the left-handed coupling of the W to the a-type parton, £ is the convolution variable 
of Eq. ([17]), mi is the initial parton mass, m-i is the heavy quark mass, and \ an d A are given by 



(Q 2 - mf + ml) + A 
~2Q 2 

A = A[-Q 2 7 mlm 2 2 ] (24) 



where rj (Eq. ([57])) is the target-mass corrected Bjorken x, and A is the triangle function (Eq. (45)) 
both defined in Appendix X[ 



Substituting in Eq. (pi), we obtain simple but non-trivial formulas for the hadron helicity struc- 
ture functions. The delta function in Eq. (22) fixes the momentum fraction variable £ = x- Since 
u>o =/= 0, we see explicitly that the longitudinal structure function cannot be neglected even to leading 
order. It is proportional to the quark masses when they are non- vanishing; thus, the Callan-Gross 
relation does not apply in its original form. 
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For charm-production, the initial parton is either a d or s quark; both can be treated as massless. 
In the limit mi — > 0, one obtains 



10+ 
10- 



= 

= g 2 La 2S^/ X -l) 



(25) 
(26) 



(27) 



and x = ^(1 + m 2/Q 2 )- Thus, the helicity structure functions assume the following simple form: 



F+ = 

F_ = g 2 La 2q a N ( X ) 

2 

^0 = gl a ^2q%(x) 



(28) 
(29) 

(30) 



where an implicit sum over contributing parton species a is implied. By applying the general 
expression of Eq. (|7]) , one obtains 



dxdy ~ G w g Le g La 2q N { X ) n 



1 + cosh ip 



m\ sinh if) 
2Q2 2 



(31) 



where tp is defined by Eq. (|63[), gn = 1 and g^ a — cos#c(sin#c) for a = s(d), respectively. Note, 
Gw = 9 2 Bw /(Q 2 + M&) = (G F /V2)/(1 + Q 2 /M&. 

The corresponding formula for anti-quark production via lepton scattering, obtained from the 
interchange of gLa and gfj a in the expressions for lo\, yields: 



F_ = 



^0 = 9la^^T N {x) 



2Q 2 



and 



da v 
dxdy 



G 2 w g 2 L£ g 2 L a2r N (x) 



vQ 2 



1 — cosh-0 



m?, sinh -0 
2Q 1 2 



(32) 
(33) 

(34) 



(35) 



These results still retain the full kinematic target-mass dependence (c/., Appendix |A|). If one 
sets M = 0, the expressions for the cross section in Eqs. ( |3l| ) and ([35]) stay unchanged; only the 
definitions of i\> and x simplify. In particular 



X — \ V 1 



M->0 



x[l + 



(36) 



which is the "slow-rescaling" variable. 



6 Comparison with Existing Calculations 

There are a variety of "slow-rescaling" prescriptions in the literature with varying degrees of accuracy. 
Some analyses of charm production in DIS use a slow-rescaling corrected parton model prescription 
which consists of using the familiar zero-mass parton model cross-section with the substitution: 

+ (37) 

This prescription incorporates only the heavy quark mass effect for the on-mass shell kinematics — 
the delta function of Eq. ([22]) — but ignores corrections to the "body" of the partonic (hard) structure 
functions lo\ in the same equation. It is therefore inherently inconsistent. 
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An improved treatment is obtained by using the exact expression for the Born diagram with 
mi = and M = 0. The results are simple enough so that the final 7712 dependence can be rewritten 
to appear as a "slow-rescaling" corrected formula, as follows: 



da 

dxdy 



_ ^2 2 2 
— 9li9l 



2Cf 
Try 



y 



X 



y) 



?(0 



y(y - 1) + 1 (1 - y) 

x 



(38) 

By definition, this modified prescription ignores target masss effects in the parton kinematics that 
are not necessarily small compared with heavy quark effects. Eq. (38) should be compared with 
Eq. ( |3l| ) which has implicit M dependence in cosh?/;, sinh-0, and x- 

Some papers include the target mass dependence of the cross section Eq. (|l3|), i.e., the term 
—Mxy/(2Ei), so that the cross section for neutrino production reads: 



da" 
dxdy 



r 2 n 2 n 2 



2Q^_ 
Try 



y + - (l-y) 



£ (Mxy\ 
x V 2Ei ) 



q(0 



(39) 

Numerically, this term has negligible effect; the —Mxyj(2E\) term does not approximate the true 
target mass dependence, and for all practical purposes, Eq. (138) and Eq. (|39) are identical at the 
< 2% level. 



c 
o 



> 

Q 
c 

0) 

o 
a. 





Fig. 4: Percent deviation of leading-order cross section between the "slow-rescaling," Eq. (p9|), and 
and complete, Eq. (31), for E v = 80GeV, m c — 1.5GeV: (a) da / dy(v + s — > c) integrated in x 
over the range x = [0.1, 0.6]; (b) da / ' dx{y + s — > c) integrated in y over the range y = [0.1, 0.8] 



We now present numerical results comparing cross-sections calculated using the complete leading 
order formula Eq. (31 ) with that using the slow-rescaling prescription, Eq. (|39|). In Fig. ||we compare 
the y and x dependence for v + TV — > /i~ + c + X for neutrino energies ranging from 50 GeV 
to 300 Gev — a reasonable range for fixed target experiments. For simplicity, we only consider the 
dominant sub-process: W + s — > c . As anticipated, for both the x and y distributions, the deviations 
decrease with increasing neutrino energy, (hence, increasing Q 2 ) since the M 2 /Q 2 and TO2/Q 2 terms 
are decreasing. The y distribution agrees well at large y, but deviates from the complete leading-order 
result by more than 25% for small y where the effects of the charm mass threshold are significant. 
The deviation of the x distribution ranges from a few percent at small x to > 25% at large x. Thus 
the difference between the conventional slow-rescaling prescription and our approach, which is based 
on the factorization theorem, are not negligible. The main source of discrepancy arises from the 
charmed quark mass mi which is only slightly larger than the target mass M; the latter should 
not be neglected if effects due to the former are significant. In particular, the momentum fraction 
variable £ = x which enters the precise formula Eq. ( |3l| ) is approximately: 



x 



1 



Q 2 



1 



x 2 M 2 



(40) 
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when m 2 ./Q 2 and M 2 /Q 2 are small, and mi = 0. In other words, the conventional "slow-rescaling" 
variable itself needs a target-mass correction. 

7 Conclusions 

The proper treatment of the effects of heavy quarks in the theoretical predictions of the differential 
cross section for deeply inelastic scattering processes is not completely solved in perturbative QCD. 
Strictly speaking, the familiar factorization theorem applies only to one scale problems, i.e., when 
either all quark masses are negligible compared to Q 2 , or when the heavy quark mass m is of the 
same order of magnitude as Q 2 . 

The recent higher order calculations of heavy quark production which exclude massive partons 
and focus on the gluon-fusion diagrams apply only to the region in which m 2 ~ Q 2 and require a 
totally different treatment of charged and neutral current processes. 

We formulate a unified approach to both types of processes that is based on the factorization 
theorem as an ansatz. We assume that the factorization theorem holds throughout the energy range 
of interest in the simple form W — f ® u. This ansatz produces the correct results in the regimes 
Q 2 ~ m 2 and Q 2 3> m 2 , and provides a smooth interpolation in the intermediate regions. We 
are able to treat both charged and neutral current processes by endowing the parton quarks with 
a mass and by not making a priori any assumptions about the relative importance of quark and 
gluon-initiated contributions. Instead, we take advantage of precisely the techniques that yield the 
proof of the factorization theorem to ensure that the final expressions conform to expectations in 
the Q 2 ~ m 2 and Q 2 » m 2 regions. 

Working towards this goal, we have presented here the general framework. In order to illustrate 
the basics of our approach, we have presented an explicit calculation of the lowest order contribution 
to the quark structure functions. However, this contribution by itself is not sufficient for proper 
phcnomcnological analysis of DIS cross sections because of the importance of quark-gluon mixing in 
sea-quark initiated processes. 

We have compared existing phenomenological analyses based on the lowest order process W+q — > 
Q, with the unified approach which retains all masses. For charged current charm production 
experiments (W + s — > c), the final state heavy quark mass m is comparable to the target mass M; 
hence, if the m-dependence is retained, then the M-dependence must also be retained for consistency. 
The m-dependence results in the well-known "slow rescaling" adjustment of the scaling variable and 
the cross section. The target mass also adjusts the effective scaling variable, and can shift the cross 
section by up to 25% for fixed-target experiments. 

For collider experiments such as the HERA e — p facility, we would like to study charged and 
neutral current production of charm and bottom quarks. Such processes fall in the intermediate 
region where the heavy quarks are neither Q 2 ~ m 2 or Q 2 3> m 2 ; hence, we must carefully take the 
mass dependence into account. 

In the second paper of this series we shall make use of the framework developed here to present 
a full next-to-leading order analysis of both charged and neutral current cross sections for deeply 
inelastic scattering. 

A Appendix I: Kinematics 

We summarize the details about the kinematics including target and heavy quark mass effects in 
this appendix. We begin with the lab frame kinematics for the overall process, and then examine 
the class of colinear frames including the Brick Wall (BW) frame. Finally, we consider the colinear 
frame for the partons, and relate the partonic quantities (including dot products) to the hadronic 
variables^ 



9 We use the metric g = {H } when necessary, but attempt to present the results in a metric independent 

fashion. 
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A.l Overall Process 

For the physical process 

*i(*i) + iV(P) - 
the following invariant variables are standard: 

P 2 = M 2 
Q 2 = 



£2(e 2 ) + x(p x ) 



(41) 



-q 

P-q 

„2 



: E\ — E2 

Q 2 



y 



2P-q 2Mv 

P-q _ v 
P~h ~ E~ x 



(42) 



where q — t\ — £2, and E\ and E2 are the laboratory energies of the incoming and outgoing leptons 
respectively. 



B 



N 




Fig. 5: Basic process for inclusive boson B(q) nucleon N(P) scattering: N{P) + B(q) — > X(Px), 
summed over the final state, X(Px) 

The components of the relevant 4-vectors in the lab frame are: 



el = 

z-2 — 
q» = 



M, 0, 0, 0] 

Ei, 0, 0, -E x 

E 2 , ~E 2 sm9, 0, -£; 2 cos^j 

+E 2 sm9, 0, -E 1 +E 2 cost 



(43) 



where, as throughout this paper, lepton masses are neglected. 

The cross section for the deep inelastic scattering process is given by the standard form: 



da 



1 



2A(s,m 2 ei ,M 2 ) 



l Af2 l dT 



(44) 



with M being the mass of the incident hadron, m£ 1 the mass of the incident lepton, and the triangular 
function 



A(o, b,c) = \Ja 2 + b 2 + c 2 - 2(ab + bc + ca) 



(45) 
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The sum and average over spins is given by 

E = — with no = # of initial spin states = < „ ' ~ . (46) 

n t ^ 1 2 for i± y ' 

spin spin 

dT represents the final state phase space, with all unobserved degrees of freedom to be integrated 
over, 

dT = dl 2 (2tt) 4 5 4 (P + £ 1 -P X - £ 2 ) dT x (47) 
with the notation (for invariant single-particle phase space) 

S - ^ =(5^ (48) 

and dTx representing the phase space factor for the hadronic final state. With the scattering 
amplitude given by Eq. (0), one can put the various pieces together to get: 

da = 2A( ^ M2) ^ V* W% di 2 dT' (49) 

Where G; = g 2 B / (Q 2 + M B . ) , the subscripts on g B . and M B . indicate the type of exchanged vector 
boson, dT' represents unintegrated hadron degrees of freedom (such as those associated with the 
production of a heavy quark), and the lepton (hadron) tensor L^ V (W W ^) is defined in Eq. (||) 
(Eq. (0)). For convenience, W and L are defined to be dimensionless; these depart from some 
historical definitions by simple factors such as M. The factor of 4irQ 2 comes from the normalization 
of W and L. 

Suppressing dT' , one obtains: 

= y -f Gl G 2 L.W (50) 
dx dy »7r 

Note that the gauge couplings of the bosons gsi appear explicitly whereas the chiral couplings of the 
leptons {gR£, gLt} and hadrons {gnh, gLh} are kept with the currents hence reside in the respective 
tensors. 

For completeness, we record the relations between various commonly used cross-sections: 

da da o da MEi v da , 

= 2^^— = 2MEfy—— : = ^——— (51) 



dxdy dxdQ 2 iy dQ 2 dv E 2 dE 2 d cos 9 

which can be easily derived using the kinematic definitions in Eq. (0). 

A. 2 The Colinear Frames 

Since the underlying physical process is actually the scattering of a space-like vector-boson on a 
nucleon, (c/., Fig. ||): 

B(q)+N(P) -> X{P X ) (52) 

it is more natural to use frames in which the 4-vectors (q, P) define the t — z plane. For parton- 
model considerations, it is convenient to specify these vectors in a general frame of this class by their 
light-cone coordinate components (x + , x, x~ ) , with x^ — (x° ± x 3 )/^/2, as: 



P» = (P+, 0, ^ 



t = {-VP + , 0, ^PT 
where P + is arbitrary, and n is defined through the implicit equation: 

Q 2 „„ 2 



(53) 



2q-P = — -<nM z (54) 
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r] represents the generalization of the familiar Bjorken-x in the presence of target mass, and it is 
related to the latter by: 

Clearly, 77 reduces to x in the zero target mass limit, 

(56) 

whereas, the general solution to Eq. (pq) is: 



1 

V 




(57) 



We shall refer to this class of frames as the colinear frames. The laboratory frame (with the negative 
z-axis aligned along q) belongs to this class; it is obtained by setting P + — M/^/2. The "infinite 
momentum frame," often used to derive the QCD asymptotic theorems, is obtained in the limit 
P + — > 00. Another useful frame in this class, used in the helicity formulation, is discussed in the 
following. 

A. 3 The Brick Wall Frame 

The Brick Wall (BW) frame is the natural "rest-frame" of the exchanged vector boson when its 
momentum q is space-like, q 2 = —Q 2 < 0, (c/. Fig. |l|). It is also one of the colinear frames — 
corresponding to setting P + = Q/{r]y/2) in Eq. (|53|), hence obtaining q° = and q 3 = — Q. In the 
cartesian coordinate system, (x°, x , x 2 , x 3 ), we have: 

5" = Q (0, 

= 2k l A P> 11 °' + M ,:,S! 

and we refer to this frame as the standard hadron configuration, Fig. ra, with 



0, 


0, 


-1 


0, 


0, 


+/3i 


0, 


0, 


-fa 



A P = A[-Q 2 ,P 2 ,P 2 } 
Pi = Q 2 -P 2 + Pl 

ft = Q 2 + P 2 -Px (59) 




Fig. 6: (a) The standard hadron configuration in {x, z}-space. Note that the hadron momenta are 
colinear with the z-axis, and the lepton momenta define the x — z plane, (b) This frame is 
related to the standard lepton configuration (Fig. fj] below) by a space-time rotation (i.e. boost) 
in the {x, t}-plane by the angle ip. 
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In this frame, the lepton momenta are given by: 

t% = § (cosh ^ sinhV>, 0, -l' 
^ (cosh-0, sinh?/>, 0, +1 



(60) 



which can be easily envisioned as being obtained from the standard lepton configuration (c/. the 
standard hadron configuration, Eq. (|58|)), Fig. 0, 



pn _ Q 

t± — 2 

pt* _ Q 

C 2 — 2 



1, 0, 0, -1 
1, 0, 0, +1 



(61) 



by a "rotation" in the (t — x) plane (really a Lorentz boost) by the hyperbolic angle ip- This is in 
analogy to the familiar CM rotation [in the [z — x) plane] between initial and final scattering states 
in a time-like situation. This is illustrated in Fig. and Fig. ^. 





Fig. 7: (a) The standard lepton configuration in {x, z}-space. Note that the lepton momenta are 
colinear with the z-axis, and the hadron momenta define the x — z plane; (b) The same fame 
seen in {x, t}-space. 



The hyperbolic cosine can be obtained from the formula: 

2P ■ {t x + e 2 ) 



cosh'0 = 



(62) 



Evaluating the scalar productions in the laboratory frame, we relate cosh'0 to the more familiar 
variables: 

E x +E 2 _ V 2 M 2 -Q 2 + 2 V (s -AP) — ► (2-y) 

M — * 



cosh i/j 



rfM 2 + Q 2 



y 



(63) 



VQ 2 + v 2 

In developing the helicity formalism (Appendix^), we encounter the "spin-1 rotation matrix" for 
the vector boson polarization vectors under the above Lorentz boost from the configuration Eq. ( |6l| ) 
(Fig. 0) to Eq. @ (Fig. |). The 3- dimensional d-matrix is: 





l+cosh?/> 
2 


— sinlii/> 


1— cosh^ 
2 






— sin hi/) 
V2 


cosh0 


+sinhi/j 
V2 


(64) 




1 — cosh'0 
2 


+sinhi/j 
V2 


l-\-coship 
2 





It is the SO(2,l) analogue of the familiar SO(3) rotation matrix. 

A. 4 Parton Kinematics in the QCD Parton Model 

In the QCD Parton Model (c/., Fig. ||), we have an initital state parton momentum ki, whose 
light-cone components in a colinear frame are: 
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(65) 



where £ is the fractional momentum carried by the parton. The momenta involved in the "hard 
scattering" consist of 

q + ki — > k x (66) 

where the final state, represented by the total momentum k x , consists of either a on-mass-shell single 
parton (for the case of the LO calculation) or a continuum of multi-parton configurations (for the 
NLO calculations and beyond). 

For the LO calculation presented in Section ^, with k x — &2 = fci + q, we can evaluate the 
argument of the delta function which enforces the on-shcll condition for the final state heavy quark: 

k 2 -m 2 = — (67) 



where 

(Q 2 - ml + mj)± A [-Q 2 , m\ , m|] 
2Q 

and J] is defined in Eq. (p7|) . The limits on £ (see below) dictate that the only physical root is 



X± = V 7^2 ( 68 ) 



£ = X = X+ (69) 

This variable reduces to the "slow-rescaling" variable x^ + m 2 ^/ Q 2 ) in the limit mi — > and M — > 0. 
Substituting Eq. ( |6^ ) in the second factor in Eq. j6^), we obtain 

^«-*"»' " A[-fci,ljl <70) 

When the final state consists of multi-partons (for NLO and beyond), the CM energy of the 
subprocess s must be greater than a threshold §th, which is equal to either m| or 4m|, depending 
on whether a single heavy quark (charged current case) or a heavy quarks-antiquark pair(neutral 
current case) is produced. Since 

S = (k 1 +q) 2 =ml-Q 2 + 2k 1 -q=(^Q 2 + ^ m 2 ^(J--?j>s th (71) 

it is easy to see that the threshold condition imposes the constraint £ > £th on the parton momentum 
fraction variable where 

c {Q 2 -mj + s th ) + A[~Q 2 , mj, s th ] 

Ztk = V ^ (72) 

(Note that for §th = m\, £ t h = X+ = X-) C> n the other hand, the condition that Pi = P + (l — £) > 
requires £ < 1. Hence, £, which is also the integration variable for the convolution in the fundamental 
factorization theorm (Eq. (p"7j)), has the following range: 

, . f . f (Q 2 - mf + s th ) + A[-Q 2 , m 2 , s th ] 

1>£>^ = ?7 ^ (73) 

We recall that rj is the generalization of Bjorken x incorporating the target mass effect. Thus the 
lower limit for £ is modified by both target mass and heavy quark mass. This aspect of mass- 
dependence has been overlooked in existing literature. 



A. 5 Dot Productions of Lepton and Parton Momenta 

In the explicit calculation of cross-sections using the contraction of lepton and hadron tensors, (c/., 
Section || and Appendix |cj) one needs the scalar products of the lepton and hadron 4-vectors. This 
calculation is subtle because the variable £ = /P+ is invariant for boosts along the z-axis, but 
not for other boosts or rotations. 
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In the BW frame, the light-cone components of the two parton momenta are: 

: * ( I' ff > ) 



(74) 

. _Q_ / n i | n m i ^ 



Using the explicit components of the lepton momenta given in Eq. ( |60| ) , it is then straightforward 
to show 

1 £ 2 / cosh + 1 77 2 ( cosh ip — 1 



1 £ 2 / coshV' — l\ 1 T) 2 ( cosh?/' + 1 



<*■•« - 55« 1—5— rd-n— f— J (76) 

To contrast the simplicity and symmetry of this group theoretic approach with a more traditional 
"brute force" calculation in the colinear frame, we compare: 

Q 2 e 2 (a - Af 2 + M 2 77) + m 2 (V - M 2 r, 2 - Q 2 t?) 

= 2#TmV) (77) 

Q 2 g 2 ( S - Af 2 - Q 2 /t?) + mfV^g + M 2 r? - M 2 ) 

(fcl '^ 2) = W + mv) (78) 

Although it is not obvious, Eq. ( f76[) and Eq. (|76] ) are identical to Eq. (|7Sl) and Eq. (|78|); however, 
the symmetries of the problem are more apparent in Eq. ([76]) and Eq. (|76|). 

In the limit of zero masses, we have the usual relations where (fci ■ £\) —* s/2 and (fci • £2) — > u/2 
with no £ dependence. However, if we wish to obtain the correct mass dependence, we must include 
the proper £ dependence in our calculation. 

Once we have (fci • £\) and [k\ ■ £2), we can use fci + 1\ = k 2 + £2 to easily compute the other 
necessary combinations via: 

(k 2 -£ 2 ) = (ki-h) 



(k 2 -£ 1 ) = (fa -l a ) + ^ z 2 L j (79) 
B Appendix II: Structure Functions and Cross-sections 

Since the precise treatment of the mass effects is emphasized in this paper, we include here some 
details on the derivation of structure function and cross-section formulas used in the text, especially 
for the less familiar helicity vertices and structure functions. 

B.l Tensor Amplitudes and Invariant Structure Functions 

We begin by recording the expression for the lepton tensor, Eq. (||). In the limit of zero lepton mass, 
it is: 

IT = £ u(£i) T» u(£ 2 ) ■ u(£ 2 ) T' 4 u(£i) = i — Tr[£x T" £ 2 T^} (80) 

W ^cJ Tig 

spin 

where rig counts the number of incoming helicity states. Using a general V-A coupling of the form, 
Eq. (|, 

r M = 7 AI [si«(i + 7 5 ) + .gL£(i-75)] (81) 
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the result is: 



1 i- 

nt <• 



)2 n t 



t x t 2 + t 2 l\ - g" 



V Q 2 



(82) 



The independent components of the hadron tensor Wnv are expressed in terms of invariant (i.e., 
Lorentz scalar) structure functions defined as (Eq. (|To|)), 



W* = 



pn p 



M 2 



Pqp, 

2M 2 



(83) 



M 2 



W 4 



P»q v + q"P v 
2M 2 



W 5 



2M 2 



W 6 



Contracting the lepton and hadron tensors and evaluating the scalar productions of the 4- vectors 
in the laboratory frame (c/., Eq. (Eq)), one obtains: 



W-L = 



\6ExE2 

n e Q 2 



9 2 +t 



2 sin" - Wi + cos 2 - W 2 



E-\ -h En o 9 
T Sin 2 W * 



(84) 



The structure functions {W4, W5, Wq] do not appear on the right-hand side of this equation because 
the dot product of q^ with the lepton tensor gives rise to a factor proportional to some com- 
binations of the lepton masses which is neglect here. Eq. (|8j|), in conjunction with Eqs. (|50|)-(51 ), 
form the bases for the derivation of the cross-section formula Eq. ( pd| ) in Section |3[ 



B.2 Helicity Vertices and Structure Functions 

We now turn to the calculation of helicity amplitudes, vertices, and structure functions. We use the 
helicity labels Ai_2 for the leptons; g\. 2 for the hadrons, and {m, n} for the bosons. Lower indices 
are for incoming particles; and upper indices are for outgoing particles. The scattering ^giplitude 
for the basic process, Eq. (|l|), can be written in the factorized form in the helicity basis: 

Mtll = J:? m {Q\t-P) ^mi" ^ 2 "(Q 2 ) (85) 

where c? 1 (V') m n is a spin-1 50(2, 1) "rotation matrix" in the Brick- Wall frame of the process cor- 
responding to q^ : ( 0, 0, 0, — Q), (c/.,Eq. (p4|)). The scalar lepton helicity vertex function 
is: 

3xl n iQ) = e^(£ 2 ,X 2 \r\£ 1 ,X 1 ) = u X2 {i 2 ) e"* T u Xl (h) (86) 

and the corresponding hadron vertex function is: 

j;Z(Q 2 ,q-P) = (Px,<7 2 |jt|P, CTl )C (87) 

Much of the simplicity of the helicity approach results from the fact that the lepton vertex function 
is extremely simple in the limit of zero lepton masses. For left-handed (right-handed) coupling, there 
is only one non-vanishing vertex function for which all three particles are left-handed (right-handed); 
it is simply given by: 

3l l {Q)=3Z\ / /1~ 1 {Q) = (88) 

(Likewise, j RR (Q) = —\/8Q 2 in the case of right-handed coupling). Thus, upon squaring the 
scattering amplitude, Eq. (jS5|), one obtains: 

53 |M 2 | oc d 1 ^)" 1 ™ dH-VO"-! W m n (89) 

spin 
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where W m n is the helicity forward Compton scattering amplitude for initial state vector boson 
polarization n and final state polarization m: 



W m n = e^(P,q)W v (P i q)^(P,q) 



(90) 



For totally inclusive process, this amplitude must be diagonal in (ro, n) due to angular momentum 
conservation;^ hence, the right-hand side becomes d 1 (ip)~ 1 m d}(ip)' m -i F m where the diagonal 
helicity amplitude W m m is identified with the helicity structure function F m , c/., Eq. (Jllf). 

Using these results for the squared amplitude, \M 2 \, keeping all factors, and making use of the 
explicit form of the d- matrix, Eq. (|64|), we obtain L ■ W, which appears in the cross-section formula 
Eq. ©: 



W-L = — < 

ni 



9 m 



F + 



1 + cosh ip 
2 

1 — cosh ip 



F 
F 



— sinh ip 
+ sinh ip 



F- 



1 — cosh ip 
2 

1 + cosh ip 



(91) 



This leads to the general formula, Eq. (0), for the cross-section given in Section^. 

B.3 Relations between Invariant and Helicity Structure Functions 

To derive the relations between the invariant and helicity structure functions, we first examine the 
polarization vectors for a vector boson with momentum q in the helicity basis. With respect to an 
arbitrary reference momentum p, the "longitudinal" polarization vector is: 



(-q 2 )p' l + (p-q) 

V(-q 2 )[(p-q) 2 -q 2 p 2 } 



(92) 



with — q 2 = Q 2 > for space-like q^. It is also useful to define the "scalar" polarization: 



<£(?,«) = 



(93) 



In a collinear frame where the z-component of q^ 1 is positive, the transverse polarization vectors are 
given by: 



4 (p. i) = ^(o»Ti,-*,o) 



(94) 



For the z-component of q M negative, we rotate the above about the y-axis by w. These polarization 
vectors depend on the reference vector p M only to the extent that it defines the t — z plane in 
conjunction with q^ . For the transverse polarization vectors, this is obvious. For the longitudinal 
vector, €q (p, q), this follows from the fact that it is merely the unit vector in the t—z plane orthogonal 
to q 1 *. The reference vector p^ is used only to define this plane and to provide the non- vanishing 
perpendicular component for projecting onto Cq. The two distinct reference vectors in the plane, 
such as P^ (the target momentum) and (the initial state parton momentum) used in the text, 
define the same set of polarization vectors for the vector boson. As discussed in Section ||, this is 
the key point which leads to the simple factorization formula for the helicity structure functions in 
the QCD Parton framework. 

To project out the transverse helicity amplitudes, the following representations are useful: 



4(p,<z) 4*(p,?) -et(p,q) e»_*(p,q) = 



q 2 p 2 



e»(p,q) e» + *(p,q) +e»(p,q) e»_*(p,q) 



-<r + <$(p,q) e»o*( P ,q) 



to rrip/Q 2 , we only concern ourselves with {W^ + j W°o, W —}■ 



$(p,q) e»*(p,q) (95) 

i W q o, W° q }. Since the coefficients of these terms are proportional 
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The second relation is simply completeness. 

Applying the above polarization vectors to the definition of the hclicity structure functions, 
Eq. (§), 

F x = e£*(P, q) Wu(P, q) *\ (P, ?) (no sum over A) (96) 

and using the representation of W ,J ' l ,(P, q) in terms of the invariant structure functions, Eq. (|8 
we obtain: 



F- = Wx + ^^1 + ^W 3 (97) 
Fo = -Wi + + W 2 

The complete transformation matrix to convert hadron helicity amplitudes to invariant amplitudes 
{W\ = / ® uj\ = t \Wi) is given in Table |[ The coefficients for the inverse transformation, (t^ 1 ) f, 
are given in Table || 

B.4 Relations Between Hadron and Parton Tensors 

As discussed in Section |], the k\ 4- vector is not simply proportional to P, but in general contains a 
mixture of P and q given by: 

Q 2 e + mW 



£(q 2 + mV) 

?7(m 2 - M 2 g 2 ) 
£(Q 2 + MV) 



Note that this mixing depends on both M and mi. The result is that the hadron tensors and the 
parton tensors are mixed. Specifically, 

W l = c{ f® uj 3 (99) 

where the c\ coefficients are given in Table [s]. The coefficients for the inverse transformation, (c _1 ) j , 
are given in Table |[ 

This is in contrast to the corresponding result for the hadron helicity amplitudes where there is 
no mixing: 

Fx = W xx = f®iu x (100) 
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Tab. 3: Transformation matrix to convert hadron helicity amplitudes to invariant amplitudes: 
W\\ — f (g) lo\ — t\Wi. Note, we use the short hand notation F\ = W\\. We have 
defined p 2 = 1 + Q 2 jv 2 , and we have p — > 1 in the DIS limit. 
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Tab. 4: Transformation matrix to convert hadron invariant amplitudes to helicity amplitudes: W, = 
(i -1 )^ Wa- Note, we use the short hand notation Fx = W\\. We have also used F\ =W\, 
F 2 = (u/M)W 2 , and F 3 = (v/M)W 3 . We have defined p 2 = 1 + Q 2 jv 2 , and we have p -> 1 
in the DIS limit. Note that asM^O, {W 4 , W 5 , W 6 } decouple from F , 
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Fi = W 1 
F 2 = {v/M)W 2 
F 3 = {u/M)W 3 

W 4 

W 5 

W 6 



1 









U>2 


U) 3 


0J4. 
















c 2 

2x 














(p 

2x 








C, 2 M 2 

Q 2 





M 2 

Q 2 


C,M 2 

Q 2 


Q 2 








CpM 2 

Q' 2 















^6 









CpM 2 



Tab. 5: Transformation matrix to convert parton invariant amplitudes to hadron invariant ampli- 
tudes: Wi = c\ f <E> ujj. Note that as M — > 0, {W4, VT5, W^} decouple from {uJi}. 



F 1 = W 1 F 2 = {v/M)W 2 F 3 = {v/M)W 3 W 4 W 5 
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Tab. 6: Transformation matrix to convert hadron invariant amplitudes to parton invariant ampli- 
tudes: ujj = (c^ 1 )) f®Wi. 



C Leading Order Calculation with Masses 

We present the details of the leading order calculation with the full mass dependence both as an 
illustration of general points made in the text of the paper, and as a concrete example to check the 
self-consistency of the tensor and helicity formalisms developed in the text. Although the calculation 
is straightforward, the results with the full mass dependence do not exist in the literature, and have 
not being used in the analysis of experimental data — as emphasized in this paper. 

The parton structure function tensor w/f", representing the vector boson (B) and parton (a) 
forward Compton scattering amplitude, is entirely analogous to W^ V N — replacing the hadron target 
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N by the parton target a. The leading order diagram, Fig. ^, gives rise to: 



An 



{2n)5+{kl - m 2 2 )J2 {kucri\f\k 2 ^ 2 ){k 2 ,a 2 \jt\k u a 1 ) 



(101) 



For quarks, the spin sum/average on the right-hand-side is: 

ilirKfr + mi)T»fa + m 2 )T v *] = Ag 2 Ra {-gT(ki • fe) + + k%h v x + ie^ p °k lp k 2a } 

+ *gl a {~9^(ki • hi) + Kk v 2 + k»k\ - ie^k lp k 2a } 



+ 4(9Ra9L a + 9L a 9R a ) {+5"* , (mim 2 )} 



(102) 



where {gn a , gL a } are the couplings of the a-type parton to the boson, and the on-mass-shell delta 
function is given by Eq. (|67| ) . 

This expression for 10^ can be used in two ways: (i) it can be substituted into the general 
factorization theorem formula, Eq. (|l7|), and then contracted with L p to yield leading order cross- 
sections directly, c/.Eq. ( [so]); or (ii) it can be used to calculate the helicity structure functions 
through Eq. ( po|) and Eq. (21) before substituting into the general cross-section formula Eq. (H). 
We shall do both, and demonstrate the consistency of the two approaches. Although at leading 
order, these two methods are comparable in the ease of use, the helicity approach provides a more 
efficient way of calculating higher orders. It also provides additional insight on the stucture of the 
physical amplitudes, as we will discuss. 

We begin with the helicity approach using 



1 

An 



2n6+(k 2 -m 2 ) 



E 



J* 



V 



no sum on A 



(103) 



and Eqs. (101)— (102) above for uj^ v (k,q), the helicity structure functions at the parton level can be 
evaluated. We obtairJPI 



-x = <M (3k n ™ + 2 9^9L a n* L + g 2 La 



(104) 



where the superscripts (R, L) refer to right-handed and left-handed chiral couplings at the hadron 
vertices, and the fi's are given in Table 0. 

The partonic helicity structure functions {u>x} exhibit many physically interesting features which 
are obscured in the conventional Dirac trace method. For example, there are obvious symmetries 
under gn a «-> g^ a when the vector boson helicity is nipped. Additionally, there is a clear order of 
magnitude separation of the amplitudes when m 2 2 /Q 2 become small (high energy limit) — all the 
longitudinal structure functions, as well as the mixed chirality ones, become of (D(m\ 2 /Q 2 )- 

Because of the direct relationship between the hadronic helicity structure functions {Fx} to the 
partonic helicity structure functions-j^}, the {Fx} functions are essentially given by the expressions 
above multiplied by the relevant parton distribution functions evaluated at £ = \ (due to the delta 
function in Eq. ( |104| )). Substituting these expressions in the general formula for L - W, Eq. (|9l|), wc 
obtain: 



L ■ W = q(0 



9 m 



9l, 



1 + cosh -0 



1 — cosh tp 



■ w 



UJQ 



sinh ip 
- sinh ip 



+ UJ- 



Note that we have used Q K = Q to simplify Eq. (104), and u> is symmetric under Q 



1 — cosh 1p 



1 + cosh lj) 



(105) 
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Tab. 7: The helicity amplitudes for the leading-order process t\ + fci(mi) — > ^2 + k 2 (m 2 ), with 
A = A[-Q 2 ,m 2 ,m|]. 



with wo, w_} given by Eq. (104). The corresponding results for the anti-quark process is 
obtained by the substitution gR a <-> <7x, o . 

Alternately, we can compute this in the tensor representation by contracting ujtf with LJ^, Eq. (|8z 
to obtain: 



L ■ u> = 



1 2 6 



(M 



n e Q 2 A[-Q 2 ,m 2 ,m 2 ] 



+ {9 2 R jli + 9lA) {ki ■ e 2 )(k 2 ■ h) 
- 9R a 9L a (9 2 Re + 9li) (mim 2 )(£i ■ £ 2 ) 



(106) 



Applying the convolution integral and inserting the scalar products between lepton and quark mo- 
menta derived in Appendix A. 5 into Eq. (106) leads to: 



L ■ W = 



q{x) 



1 2 b 



n e Q 2 A[-Q 2 ,m 2 ,m 2 ] 



{9rA + 9l a 9li) (Q 2 X 2 d- + m 2 r] 2 d + 
(gRjle + 9 2 L a 9 2 Ri) (Q 2 X 2 d+ + m 2 r/ 2 cL 
9R a gL a (g 2 R i + 9li) (mim 2 )Q 2 /2 



xvQ 2 )(Q 2 x 2 d, 
xvQ 2 )(Q 2 x 2 d- 



+ m 2 l rj 2 d-)/{2 2 rj 2 x 2 ) ' 
+ m 2 1 rj 2 d + )/(2 2 -q 2 x 2 ) 



(107) 



where d± — (cosh?/' ± l)/2 are elements of the d 1 ^) matrix. A special case of these results - charm 
production in neutrino scattering - is discussed in Secti on |q. 

Although it is far from obvious, Eq. (105) and Eq. ( |107[ ) are in fact identical (as some tedious 
algebra will prove). The difference in appearance is simply that the helicity approach exploits the 
symmetries of the problem; hence, these symmetries are manifest in the final representation of the 
cross section, Eq. (105). 



26 



Acknowledgement 

The authors would like to thank Andrew Bazarko, Raymond Brock, John Collins, Sanjib Mishra, 
Michael Shaevitz, and Davison Soper for useful discussions. This work is partially supported by 
the National Science Foundation under Grant No.PHY89-05161, the U.S. Department of Energy 
Contract No. DE-FG05-92ER-40722, and by the Texas National Research Laboratory Commission. 
M.A. and F.O. also thank the Lightner-Sams Foundation for support. F.O. is supported in part by 
an SSC Fellowship. 

References 

[1] Wu-Ki Tung, et al, Proceedings of the Summer Study on High Energy Physics in the 1990s, 
June 27- July 15, 1988, Snowmass, Colo.; High Energy Physics in the 1990s, (ed., S. Jensen), 
p. 305 (1989). 

[2] J. Botts, et al, Phys. Lett. B304, 159 (1993). 

[3] J.G. Morfin and Wu-Ki Tung, Z. Phys. C52, 13 (1991). 

[4] Joseph F. Owens and Wu-Ki Tung, Ann. Rev. Nucl. Part. Sci. (1992). 

[5] K. Lang, et al., Z. Phys. C33, 483 (1987); 

B. A. Schumm,e£ al, Phys. Rev. Lett., 60, 1618 (1988); 

S.R. Mishra, et al, Moriond 1989, Leptonic:3 (QCD161:R4:1989:V.l); 
Edward Oltman, Ph.D. Thesis, Columbia University report Nevis 270, 1989; 
Coustantinos Foudas, Ph.D. Thesis, Columbia University report Nevis 272, 1989; 

C. Foudas, ,et al, Phys. Rev. Lett, 64, 1207 (1990); 
P.G. Reutens, et al, Z. Phys. C45, 539 (1990); 

Michael H. Shaevitz, Nucl. Phys. B (Proc. Suppl.) 19, 270 (1991). 

S.R. Mishra, et al, Columbia Univ. preprint NEVIS-1466, (1992). Invited talk given at Lepton 
Photon '91 Conf., Geneva, Switzerland, 1991. 

[6] D. Bogert, et al, Phys. Rev. Lett, 55, 1969 (1985); 
Phys. Rev., D43, 2778 (1991); 

T. S. Mattison, et al, Phys. Rev., D42, 1311 (1990). 

[7] H. Abramowicz, et al, Phys. Rev. Lett., 57, 298 (1986); 
Z. Phys. C28, 51 (1985); 
J.P. Berge, et al, Z. Phys. C49, 187 (1991). 

[8] U. Amaldi, et al, Phys. Rev. D36, 1385 (1987); 

See also J. Feltesse, Proceedings of the 1989 International Symposium on Lepton and Photon 
Interactions at High Energies, Stanford, August 1989, Ed. M. Riordan, p. 13, (1990). 

[9] J. Ellis and G.L. Fogli, Phys. Lett. B232, 139 (1989); 
G.L. Fogli and D. Haidt, Z. Phys. C40, 379 (1988). 

[10] Francis Halzen and Duncan A. Morris, Phys. Lett. B237, 107 (1990). 

[11] Sec, for instance, W. Furmanski and R. Pctronzio, Z. Phys. C, 11, 293 (1982). 

[12] J.J. van dcr Bij and G.J. van Oldenborgh, Z. Phys. C51, 477 (1991). 

[13] B. Lampc, Z. Phys. C34, 523 (1987); 

G. Kramer and B. Lampc, Z. Phys. C54, 139 (1992). 



27 



[14] E. Laenen, S. Ricmcrsma, J. Smith, and W.L. van Neerven, SUNY-Stony Brook preprint ITP- 
SB-92-35, (1992); 

E. Lacncn, S. Riemersma, J. Smith J. Smith, and W.L. van Neerven, SUNY-Stony Brook 
preprint ITP-SB-92-09, (1992); 

J. Smith and W.L. van Neerven, Nucl. Phys. B374, 36 (1992); 

S. Ricmcrsma, J. Smith, and W.L. van Neerven, Phys. Lett. B282, 171 (1992); 

E. Laenen, J. Smith, and W.L. van Neerven, Nucl. Phys. B369, 543 (1992); 

W. Beenakker, W.L. van Neerven, R. Meng, G.A. Schuler, and J. Smith, Nucl. Phys. B351, 
507 (1991); 

W. Beenakker, H. Kuijf, W.L. van Neerven, and J. Smith, Phys. Rev. D40, 54 (1989). 

[15] J.V. Allaby, et al, Z. Phys. C36, 611 (1987); 
Phys. Lett. 197B, 281 (1987); 
Phys. Lett. B213, 554 (1988). 

[16] M. Arneodo, et al, Z. Phys. C35, 1 (1987); 
Nucl. Phys. B333, 1 (1990); 
Nucl. Phys. B321, 541 (1989); 
J.J. Aubert,ei al, Phys. Lett. 94B, 101 (1980); 
Phys. Lett. 167B, 127 (1986); 
Nucl. Phys. B293, 740 (1987); 

[17] R. M. Barnett, Phys. Rev. Lett., 36 1163 (1976). 
R. M. Barnett, Phys. Rev. D14, 70 (1976); 
R.J.N. Phillips, Nucl Phys. B212, 109 (1983). 

[18] Thomas Gottschalk, Phys. Rev. D23, 56 (1981). 

[19] R. Brock, Talk delivered at the New Directions in Neutrino Physics at Fermilab workshop, 
Batavia, Illinois, September 1988. 

R. Brock, C.N. Brown, H.E. Montgomery, M.D. Corcoran, Brcckenridgc Workshop, p. 358, 
(1989), (QCD161:W638:1989). 

[20] M.A.G. Aivazis, Fredrick I. Olness, and Wu-Ki Tung, Phys .Rev. Lett. 65, 2339 (1990). 

[21] A. Ali, F. Barreiro, J.F. de Troconiz, G.A. Schuler, and J.J. van der Bij, Workshop on Large 
Hadron Colliders, Aachen, Germany, Oct 4-9, 1990. CERN report 90-10, p.917 (1990); 
A. Ali, G. Ingelman, G.A. Schuler, F. Barreiro, M.A. Garcia, J.F. de Troconiz, R.A. Eichler, 
and Z. Kunszt, DESY HERA Workshop (QCD161:H15:1987) p. 395, (1987). 

[22] G. Ingelman and G.A. Schuler, Z. Phys. C40, 299 (1988); 
Gerhard A. Schuler, Nucl Phys. B299, 21 (1988). 

[23] J. Collins, D. Sopcr and G. Sterman, Nucl Phys., B263, 37 (1986) 

[24] J. Collins, F. Wilczck, and A. Zee, Phys. Rev. D18, 242 (1978). 

[25] Fredrick I. Olness and Wu-Ki Tung, Warsaw Symp. (QCD161:W3:1989) p. 375 (1989); 

F. I. Olness and Wu-Ki Tung, Nucl. Phys. B308, 813 (1988). 

[26] M.A.G. Aivazis, John C. Collins, Fredrick I. Olness, and Wu-Ki Tung, "Leptoproduction of 
Heavy Quarks II - A Unified QCD Formulation of Charged and Neutral Current Processes 
from Fixed-target to Collider Energies", Preprint MSU-HEP 93/17 (to be published). 

[27] T. P. Cheung and Wu-Ki Tung, Phys. Rev. D3, 733 (1971); 
P. H. Frampton and Wu-Ki Tung, Phys. Rev. D3, 1114 (1971). 

[28] Fredrick Olness and Wu-Ki Tung, Phys. Rev. D35, 833 (1987). 

[29] John C. Collins, Davison E. Soper, George Sterman, Perturbative QCD, A.H. Mueller, ed., 
WorldScientific Publ., 1989. 



28 



